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A  STUDY  OF  SCATTERING  FROM  E-PLANE 
CIRCUIT  ELEMENTS 


This  study  contains  two  parts.  In  part  one,  a  new  analytical  technique  is 
introduced  to  characterize  the  scattering  phenomena  of  a  number  of  planar  E-plane 
obstacles.  This  technique  is  based  on  the  spectral  domain  method  combined  with 
residue  calculus.  This  is  the  first  time  the  spectral  domain  method  has  been  used  to 
solve  a  waveguide  excitation  problem.  It  provides  some  attractive  features.  For 
instance,  it  can  handle  a  wide  non-touching  E-plane  fin  on  a  dielectric  substrate;  the 
method  deals  with  inhomogeneous  algebraic  equations  instead  of  integral  equations; 
and  also  the  spectral  domain  Green's  functions  have  simple  closed-form 
expressions.  The  calculated  results  for  a  simplified  case  arc  compared  with  existing 
data,  and  are  in  good  agreement.  Several  curves,  useful  for  E-plane  configurations, 
are  included. 

Part  two  presents  an  analysis  and  design  of  an  evanescent  mode  bandpass 
waveguide  filter  with  non-touching  E-plane  fins.  The  theoretical  analysis  is  based 
on  the  generalized  scattering  matrix  technique  in  conjunction  with  the  spectral 
domain  approach,  which  is  described  in  part  one,  and  mode  matching  method.  The 
method  used  in  this  study  takes  into  account  the  dominant  as  well  as  the  higher  order 
mode  effects.  The  measured  filter  responses  in  Ka  band  (26.5-40.0  GHz)  are  in 
good  agreement  with  those  obtained  by  this  analysis. 
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PART  ONE 


SPECTRAL  DOMAIN  ANALYSIS  OF  SCATTERING 
FROM  E-PLANE  CIRCUIT  ELEMENTS 


l 


I.  INTRODUCTION 


Recently,  finlines  [1,2]  and  other  E-plane  structures  [3,4,5]  have  found 
wide  applications  in  millimeter-wave  integrated  circuits.  Their  advantages  over 
conventional  microstrip  at  millimeter  wavelengths  include:  reduced  radiation,  less 
stringent  tolerances,  compatibility  with  hybrid  devices,  and  ease  of  transition  to 
standard  waveguide.  A  number  of  passive,  active  and  non-reciprocal  components 
have  been  developed  with  the  E-plane  technique.  One  of  the  key  elements  for 
passive  E-plane  components  is  the  E-plane  strip  that  is  inserted  in  the  middle  of  a 
waveguide  parallel  to  the  E  plane,  as  shown  in  Fig.  1-1.  A  comprehensive  design 
process  of  E-plane  bandpass  filters  has  been  reported  [6].  The  E-plane  fin  described 
in  Fig.  1-1  is  an  inductive  element.  Its  equivalent  circuit  is  shown  in  Fig.  I-l(c). 
The  analysis  of  the  E-plane  fin  connecting  the  top  and  bottom  walls  is  relatively 
straightforward  [6],  because  the  problem  is  two-dimensional.  The  non-touching 
E-plane  fin  described  in  Fig.  1-2  is  the  more  general  form.  It  contains  both  inductive 
and  capacitive  elements,  as  shown  in  its  equivalent  circuit.  It  is  more  flexible  for 
design.  On  the  other  hand,  no  extensive  and  accurate  characterizations  of 
non-touched  E-plane  fins  seem  to  exist.  A  method  based  on  a  variational  technique 
has  been  introduced  for  a  special  case  where  there  is  no  dielectric  substrate  inserted 
in  the  waveguide  [7].  The  method  in  [7]  is  useful  for  a  narrow  strip,  because  only 
one  current  component  along  the  E-plane  direction  is  used  and  the  assumed  current 
distribution  is  constant  in  the  axial  (Z)  direction. 

This  study  introduces  a  new  analytical  technique  to  characterize  the 
scattering  phenomena  of  a  number  of  planar  E-plane  obstacles.  For  instance,  it  can 
handle  a  wide  non-touching  E-plane  fin  on  a  dielectric  substrate.  Unlike  the  method 


based  on  the  variational  technique,  scattering  coefficients  of  the  dominant,  as  well  as 
higher-order,  modes  can  be  derived.  The  incident  mode  can  be  either  dominant  or  a 
higher-order. 

The  method  in  this  study  is  an  extension  of  the  spectral-domain  method 
commonly  used  for  characterizations  of  eigenmodes  in  a  transmission  line  on  a 
multi-layer  structure.  It  is  extended  to  the  excitation  problem  and  hence  provides  a 
set  of  algebraic  equations  corresponding  to  coupled  integral  equations  that  would  be 
derived  in  the  space  domain.  Compared  to  the  integral  equation  method,  the  new 
technique  has  a  number  of  advantages.  For  instance,  algebraic  equations  are  easier 
to  handle  numerically.  Also  the  spectral-domain  Green's  functions  have  simple 
closed-form  expressions.  Compared  to  the  variational  method  [7],  the  present 
method  is  not  only  more  versatile  but  also  is  attractive  from  a  computational  point  of 
view.  In  the  new  method,  it  is  necessary  to  calculate  the  eigenvalue  of  only  the 
particular  scattered  mode  of  interest  The  variational  method  requires  evaluations  of 
all  eigenvalues.  Furthermore,  the  method  in  [7]  assumes  that  only  the  TE-  modes 
are  scattered.  By  the  nature  of  the  formulation,  the  present  method  contains  both  the 
TE-  and  TM-  modes  in  its  formulation. 

The  calculated  results  for  a  special  case  (Ej=e  2  =£  3=  1.  See  Fig.I-2  (b)) 

are  compared  with  experimental  and  computed  data  in  [7]  to  check  accuracy. 
Several  useful  curves  of  normalized  input  admittance  and  equivalent  circuit  element 
values  are  presented  for  a  number  of  different  parameters  of  the  structure.  Those 
data  are  applicable  to  the  design  of  bandpass  filters,  diode  mounts,  and  tuning 
elements  of  the  form  described  by  Konishi  [3,4,5]. 


n. 


FORMULATION 


In  this  study  only  the  unilateral  E-plane  fin  is  treated.  Application  to  a 
bilateral  configuration  is  straightforward.  In  reference  to  Fig.  1-2,  the  strip  is 
assumed  to  be  perfectly  conducting  and  infinitesimally  thin.  The  first  step  of 
analysis  is  to  find  the  scattering  parameters  at  the  edge  of  the  non-touching  E-plane 
fin.  Then  the  normalized  admittance  and  other  equivalent  circuit  element  values  can 
be  obtained  from  the  scattering  parameters.  The  basic  idea  used  in  the  solution  of 
this  problem  is  described  briefly  as  follows. 

For  a  given  incident  field  E1,  the  Fourier  transform  of  the  field  scattered 

by  the  strip,  say  Ey,  can  be  obtained  by  using  the  spectral  domain  technique.  On  the 

other  hand,  the  scattered  field  E  can  also  be  written  as  a  sum  of  the  eigenfunctions 

with  unknown  coefficients.  Taking  the  Fourier  transform  of  the  sum  and  applying 
the  orthogonality  of  the  eigenfuctions,  one  determines  the  unknown  coefficients, 
which  are  related  to  the  scattering  parameters.  The  non-touchting  fin  can  be 
expressed  by  its  equivalent  circuit.  When  the  waveguide  is  terminated  with  a 
matched  load,  the  values  of  the  equivalent  circuit  elements  are  calculated  from  the 
real  and  imaginary  parts  of  the  normalized  input  admittance.  The  procedure  is  given 
in  detail  below. 

A.  MODAL  EXPANSION  OF  SCATTERED  FIELDS 

In  this  study  the  harmonic  time  dependence  is  the  form  in  e  J051,  but  this 
term  is  suppressed  throughout  the  analysis  for  convenience.  Consider  an  incident 


electric  field  E‘  which  comes  from  left  side  of  the  waveguide  as  indicated  in  Fig.I-2, 


the  modal  expansion  for  the  Ey  component  of  the  scattered  field  is  given  by  : 
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Here  eQ  and  |iQ  are  the  permittivity  and  permeability  of  free  space,  respectively;  £j 


is  the  relative  dielectric  constant  in  the  region  i;  an  is  the  nth  eigenmode 

corresponding  to  the  Y  direction  in  the  waveguide;  Yjm  and  Yifn  are  the  mth 

eigenmodes  of  LSE-  and  LSM-modes  in  the  partially  filled  waveguide  in  region  i, 
and  can  be  obtained  by  solving  the  following  eigenvalue  equations  [8],  together  with 
equations  ( 1-5)  and  ( 1-6) . 


a 


tanh  (Yjh,) 


^  Y^tanh  (Y;h2) 
YBnh(vh3) 


=  0 


tanh  (Y3h3) 


+  Y2tanh(Y2h2) 


(1-8) 


8 


•w 

S3 


YjtanMYjh,)  i 


Y3tanh(y3h3)  Y2«anh(y2h2) 


£2  1  Y3tanh(Y2h2)tanh(y3h3) 


=  0 


(1-9) 


h  1ft 

where  pmn  and  P'^  are  the  propagation  constants  of  the  ninth  LSE-  and 
LSM-mode,  respectivelty.  An  equation  similar  to  (  1-1  )  can  be  written  for  Ez 
component .  If  the  induced  current  components  Jy(y,z)  and  Jz(y,z)  are  provided  the 
scattered  fields  Ey  and  Ez  can  be  expressed  as  follows  : 


fit 
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where  Gyy  etc.  are  the  Green's  functions  in  the  space  domain.  One  way  to  find  Jy 

and  Jz  is  the  application  of  the  integral  equations  which  require  that  the  total 
tangential  electric  field  components  to  be  zero  on  the  strip. 

Ey(x,y,z)  +  E?(x,y,z)  -0  (1-12) 

(x  -  hx  +  h2,  y,  z  on  strip). 
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The  integral  equations  are 
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If  Eqs.  (1-14)  and  (1-15)  are  solved,  Jy  and  Jz  can  be  obtained.  These  Jy  and  Jz  are 

then  substituted  into  Eqs.  (1-10)  and  (I- 11)  so  that  Ey  and  Ez  are  available 

everywhere.  If  the  scattered  field  coefficient  of  a  particular  mode  is  needed,  the  Ey 
$ 

and  Ez  can  be  used  in  ( 1-1 ).  Each  coefficient  may  be  found  using  the  orthogonality 

of  the  expansion  functions.  Although  the  above  formulation  is  correct,  such  an 
approach  is  not  adopted  in  this  study.  Instead,  a  corresponding  procedure  is  applied 
in  the  Fourier-transformed  domain.  There  are  two  reasons  for  using  this  new 
technique. 

(1)  In  the  Fourier-transformed  (spectral)  domain,  coupled  algebraic 
equations  instead  of  coupled  integral  Eqs.  (1-14 )  and  (1-15)  are  involved. 

(2)  Derivation  of  the  Green's  functions  in  the  space-domain  is  very 
complicated.  In  the  spectral-domain,  the  Fourier-transformed  Green's  functions  are 
given  in  closed  forms. 

The  Fourier  transform  is  defined  as  following: 

F(a„,P)  =  eim’y dy f°°  f(y,z)e^dz 

J -  b  •'-oo 


The  Fourier  transform  of  Eq.  (1-1)  at  x=h j+h2  is  given  by 
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where 


_/2  n-< 
1  \l 


fl  n  +  0 
1  “  10  n-0. 


The  coefficients  C^n  and  D^,n  can  be  solved  by  the  orthogonality  of  the 
eigenfunctions,  if  the  scattering  field  E*  is  known. 


B.  SCATTERED  FIELDS  IN  SPECTRAL  DOMAIN 


Next  step  is  to  obtain  the  scattered  fields  in  the  Fourier-transformed 
domain.  To  this  end,  the  spectral-domain  technique  needs  to  be  extended  to  the 

excitation  problem.  The  Fourier  transform  of  Eqs.  (1-14)  and  (1-15)  at  x=hj+h2  is 
given  by 

£'(«.,^)+ d-19) 

where 


£/(<*«. 0)  “^(a..^)4(®«.^)  +  <^,(®..^)-4(«.,/?) 


are  the  scattered  electric  fields  in  the  spectral  domain.  Gyy,  G^,  and  Ga  can  be 

obtained  by  the  immittance  approach  [9]  ( Appendix  A).  Notice  that  the  right  sides 
of  Eqs.  (1-18)  and  (1-19)  are  not  zero.  This  is  because  the  application  of  the  Fourier 
transform  requires  use  of  fields  not  only  on  the  strip  but  also  the  fields  outside  of 

tthe  strip.  Hence,  Eqs.  (1-18)  and  (1-19)  contains  four  unknowns  Jy,  Jr  Ey  and  E^. 
In  the  process  of  solution  by  Galerkin's  method,  E^  and  E^  are  eliminated.  To  this 
end,  Jy  and  Jz  are  expand  as  follows: 


(1-23) 


where  *  indicates  the  complex  conjugate. 


For  a  given  incident  field,  Eqs.(l-24)  and  (1-25)  are  solved,  and  a^  and  bj 

can  be  obtained.  Hence,  Jy  and  Jz  are  now  known.  The  scattered  fielts  in  the 
Fourier-transformed  domain,  which  is  described  by  Eqs.  (1-20)  and  (1-21)  ,  now  are 
completely  determined.  Since  Gyy ,  Gyz ,  ,  Jy  ,and  Jz  are  all  known. 


C.  SCATTERING  PARAMETERS 


The  scattering  coefficients  C^n  and  D^n  can  now  be  obtained.  Let  us 


express  the  left-hand  side  of  Eq.  (1-16)  with  Eq.  (1-20).  Since  Jv  and  J,  are  now 

y  z 


found,  Gyy  and  Gyz  are  given  in  closed  forms:  (Appendix  A) 


G  —  Z'N,1  +  ZhNl 


(1-31) 


<v-$,-(z'-z»Kiv, 


(1-32) 


G„  -  ZV/,1  +  Z'N* 


(1-33) 


where  Ze,  Zh,  Ny,  and  Nz  are  given  in  Appendix  A.  The  left-hand  side  of  Eq.(I-17) 


is  completely  known.  Furthermore,  the  left-hand  side  contains  poles  at  p=-p  and 


P  =iPmn>  since  they  are  zeros  of  the  denominators  of  Gyy  and  Gyz.  These  values 


provide  the  eigenvalues  of  the  LSE-  and  LSM-modes,  yim  and  y'm.  The  right-hand 


16 


side  of  (1-17)  contains  LSE  poles  at  -($mn  in  the  C^n  term  and  at  +(3mn  in  the  C+ 


'mn 


mn 


'mn 


term.  It  contains  LSM  poles  at  -|5'mn  in  the  D^n  term  and  at  +(3'mn  in  the  D, 


mn 


mn 


mn 


term.  Therefore,  C^n  and  Dj^n  can  be  obtained  by  residue  calculus.  This  process 

corresponds  to  the  use  of  the  orthogonality  relationship  to  find  the  modal  coefficients 
in  the  space-domain.  The  results  are  given  by: 


C*.m±  lim 
P~  ±Pm 


^  >66asinh(YJm/ij) 


>*-0.1.2.  -/?  n"°) 

"  11  n  +  0/ 


(1-34) 
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(1-35) 


It  should  be  noted  that,  the  above  equations  are  simplified  in  the  case  of  E-plane  fins 
connecting  the  top  and  bottom  wall.  Since  there  are  no  field  variations  in  the 

y-direction  and  only  TEno  modes  are  scattered  for  a  TE,0  excitation,  we  have  only 


one  equation  (1-18).  All  the  Fourier- transformed  quantities  are  functions  of  (3  only. 


D.  NORMALIZED  INPUT  IMPEDANCE  AND  EQUIVALENT  CIRCUIT 
ELEMENTS 


The  non-touching  E-plane  fin  shown  in  Fig.  I-2(a)  may  be  represented  by 
an  equivalent  T-network,  as  shown  in  Fig.  I-2(c).  When  the  waveguide  is 
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terminated  with  a  matched  load,  as  shown  in  Fig.  I-3(a),  its  equivalent  circuit  is 
shown  in  Fig.  I- 3(b).  The  normalized  input  admittance  may  be  represented  by 


1  R  .  . 

y.  - G  +  iB  - -  (1-36) 

“  i*  -'m  R 

where  R  is  the  reflection  coefficient  for  the  dominant  mode,  and  can  be  determined 
in  the  present  method  by  means  of  Eqs.  (1-34  )  and  (1-35).  The  normalized 

reactances  X!  and  X2can  be  determined,  after  Yin  in  Eqs.  (1-34)  and  (1-35)  is 

obtained.  The  expressions  of  Xj  and  X2  are  given  by  Eqs.  (1-37)  and  (1-38). 
(APPENDIX  B) 


[( 

U£+<?£-<?J 

Bl  +  Gl-G*' 
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(1-37) 

(1-38) 


The  sign  in  Eqs.(I-37)  and  (1-38)  can  be  determined  by  using  the  Foster 

reactance  theorem  which  require  dx/do)  >  0  for  a  lossless  element.  We  can  also 
determine  the  values  of  C  and  L  in  the  equivalent  circuit  Fig.  I-3(b),  once  X2  is 
found,  under  the  assumption  that  the  variations  of  C  and  L  with  frequency  are  small. 
The  capacitance  and  inductance  can  be  obtained  by  solving  coupled  equations. 


in  RESULTS  AND  DISCUSSION 

A.  CHOICE  OF  BASIS  FUNCTION 

In  principle  any  kind  of  basis  functions  may  be  used  as  long  as  it  is 
nonzero  only  on  the  strip.  However,  due  to  the  variational  nature  of  the  spectral 
domain  method,  the  efficiency  and  accuracy  of  this  approach  depends  greatly  on  the 
choice  of  basis  functions.  For  rapid  convergence  of  the  solution,  the  functions 
should  satisfy  the  edge  condition  [10].  Also  it  is  desirable  that  the  Fourier 
transforms  of  the  basis  functions  are  available  in  a  analytical  form.  With  the  above 
consideration,  the  following  set  of  functions  are  employed. 


■MS-!)]  -MS*?)] 


where  k  and  /  are  integers,  and  i  and  j  are  given  by  a  suitable  combination  of  k 
and  /.  The  basis  function  Jyl(y,z)  (/  =1  and  k  =1  )  is  shown  in  Fig.  I-4(a);  and 

Jzi(y.z),  corresponding  to  /  =1  and  k  -2,  is  expresed  in  Fig.  I~4(b  ).  The  Fourier 
transforms  of  the  basis  functions  are  given  by  Eqs.  (1-43)  and  (1-44). 
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Jyi  (otj.p)  =  (  0.5d  e  0.5j(2i-i)W  { jQ  ,  o.5(2 ft- 1)jc  +  a,  d  I 

-  e-ie*-1*  Jc  I  a,  d  -  0.5(2ft-l)jt  I  ]  }  {0.25wtc  e  0-5j(k-l)n 

[  JD 1 0.5(k-l)jc  +  0.5  Pw  I  +  c-K2l'l)«  JQ 1 0.5  Pw-  0.5(k-l)7t  I  ]} 

(1-43) 

Jzj(<Xj,p)  =  {  0.5d  e  03j(2ft-i)n  {  Jq  ,  o.5(2l  -1)tc  +  otj  d  I 

+  JQ I  a,  d  -  0.5(2ft-l)7C  I  ]  }  (0.25w7t  e  0-5j(k-l)rt 

[  JD 1 0.5(k-l)jt  +  0.5pw  I  -  e*^2**1^  J0 1 0.5pw-  0.5(k-l)7tl  ]} 

(1-44) 

where  Jc  denodes  the  zero-order  Bessel  function  of  the  first  kind. 

B.  VALIDITY  CHECK 

To  compare  the  present  method  with  the  experimental  and  computed  data 
in  the  previous  publication  [7],  we  considered  first  special  case  e^e^e^l.  We 

assume  that  a  dominant-mode  incident  electric  field  Ey10  =  -jP10<|)1(x)e'->PloZ, 

where  <J>j(x)  is  defined  in  Eq.  (1-2),  comes  from  the  left  of  the  waveguide.  The 
numerical  result  has  been  checked  by  the  power  conservation  law  in  which  the 
equation  IRI2  +  ITI2  =  1  has  to  be  satisfied,  R  and  T  are  the  reflection  and 

transmission  coefficients  for  the  dominant-mode,  and  can  be  determined  from  Cjq 
and  Cjq,  which  are  given  by  Eq.  (1-34).  The  calculated  results  for  different 
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parameters  are  given  in  Table  1.  It  shows  that  IRI  2  +  ITI  2  is  essentially  1.  The 
convergence  test  with  different  number  of  basis  functions  has  been  performed. 
Although  accuracy  can  be  increased  with  a  larger  number  of  basis  functions, 
computation  efforts  also  increase. 

A  comparison  of  the  normalized  susceptance  versus  frequency  between 
the  numerical  results  obtained  by  present  approach  and  those  given  in  [7]  is  shown 
in  Fig.  I-5(a)  and  (b).  It  is  found  that  the  numerical  results  agree  well  with  the 
experimental  data  and  Chang's  data.  It  is  believed  that  the  present  method  is  more 
accurate  and  can  be  improved  systematically  with  the  use  of  more  basis  functions. 

C.  RESULTS  AND  DISCUSSION 

Figures  I-6(a)  and  I-6(b)  show  the  variations  in  normalized  admittance 
versus  frequency  with  different  values  of  dielectric  constant  of  the  substrate  in 
Region  2.  As  the  dielectric  constant  of  the  substrate  increases,  the  resonant 

frequency  at  which  the  Bjn  becomes  zero  decreases.  This  phenomenon  happens 

because  as  the  dielectric  constant  of  the  substrate  increases  the  wavelengths 
corresponding  to  each  mode  become  shorter.  The  resonant  frequency  and  the 
characteristics  of  input  admittance  of  the  structure  can  be  controlled  by  the  dielectric 
constant  of  the  substrate. 

Figures  1-7  and  1-8  show  the  variations  of  normalized  admittance  versus 
the  height  and  width  of  the  strip  at  different  frequencies  for  an  E-plane  fin  inserted  in 
a  Ka-band  rectangular  waveguide.  As  shown  in  Fig.  1-7,  for  a  fixed  width  of  the 
fin,  the  shorter  the  height ,  the  higher  the  resonant  frequency  is,  because  the  shorter 


fin  corresponds  to  smaller  values  of  capacitance  and  inductance  in  the  equivalent 
circuit.  For  a  fixed  height  of  the  fin  (Fig.  1-8),  the  higher  resonant  frequency  is 
related  to  a  narrower  fin,  which  has  smaller  capacitance.  One  way  to  obtain  a  given 
resonant  frequency  of  this  structure,  is  to  change  the  dimension  of  the  fin. 

Figures  1-9  and  I- 10  show  the  values  of  Xj  and  X2  versus  the  height  d 
and  the  width  w  of  the  strip  in  an  X-band  rectangular  waveguide.  It  is  noted  that  for 
a  given  frequency,  Xj  is  not  sensitive  to  d,  whereas  X2  increases  with  d.  When  w 

increases,  Xj  increases.  On  the  other  hand,  X2  decreases  with  w  for  a  higher 
frequency  and  increases  in  a  certain  region  of  w  for  lower  frequency.  (Curve  A  in 
Fig.  I- 10(b)).  Figures  1-11  and  1-12  are  the  data  corresponding  to  those  in  Figs.  1-9 
and  I- 10  except  that  the  frequencies  are  in  the  Ka-band.  Figures  1-13  and  1-14  show 

Xi  and  X2  versus  frequency  with  different  values  of  height  and  width  for  the 
Ka-band  waveguide.  In  Fig.  1-13  it  is  seen  that  for  a  narrow  strip  Xj  varies  slowly 

as  frequency  is  increased,  while  X2  varies  faster  for  a  narrower  strip  than  for  a  wider 
one.  Those  informations  are  very  useful  from  design  point  of  view.  For  instance, 
to  control  the  resonant  frequency  of  this  structure,  it  is  much  better  to  tune  the  height 

of  the  fin  instead  of  changing  its  width.  In  this  way,  X2  can  be  changed  without 

affecting  Xj  since  only  X2  is  sensitive  to  the  variation  of  the  height  of  the  fin  . 

Figure  1-15  shows  the  variation  of  normalized  capacitance  C  and 
normalized  inductance  L  versus  d.  It  is  noted  that  there  are  two  regions.  In  one  of 


them  (approximately  corresponding  to  d=w,  in  this  calculation  w=  1.0mm)  C 
increases  and  L  decreases  as  frequency  increases.  In  another  region  C  decreases  and 
L  increases  as  frequency  becomes  higher.  It  is  conjectured  that  this  phenomenon  is 
related  to  the  field  distributions  and  to  the  use  of  the  equivalent  circuit  chosen  here. 
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Fig.  1-12  Normalized  reactance  of  equivalent  circuit  of  an  E-plane  strip  in  a 
Ka-band  waveguide  versus  strip  width  for  e(=e^=l,  £2=2.2, 
h|=3.43mm,  h2=0. 254mm,  d=2.45,  (a)  Normalized  X]  versus  strip 
width,  (b)  Normalized  X,  versus  strip  width. 


Fig.  1-15  Normalized  equivalent  circuit  element  values  C  and  L  of  a  narrow  E-plane 
strip  as  a  function  of  the  strip  height.  e,=£j=l,  t^=2.2,  hj=3.43mm, 
h2=0.254mm,  w=lmm,  (a)  Normalized  C  versus  strip  height. 

(b)  Normalized  L  versus  strip  height. 


IV  CONCLUSIONS 


A  new  analytical  technique  has  been  developed,  that  can  be  used  for 
characterizations  of  the  scattering  phenomena  of  planar  E-plane  obstacles.  The 
numerical  results  for  a  special  case  agree  well  with  experimental  and  published  data. 
The  curves  of  normalized  input  admittance  and  equivalent  circuit  element  values  are 
presented  for  a  number  of  different  parameters  of  this  structure.  This  technique  is 
believed  to  bo  useful  in  the  design  of  microwave  filters  and  other  planar  circuit 
components. 
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PART  TWO 


ANALYSIS  AND  DESIGN  OF  EVANESCENT  MODE 
WAVEGUIDE  BANDPASS  FILTER  WITH 
NON-TOUCHING  E-PLANE  FINS 


I 


INTRODUCTION 


In  the  late  1950’s  Jaynes  [1]  and  Edson  [2]  proposed  that  resonators  built  in 
cutoff  waveguide  may  be  used  in  filter  design.  These  kinds  of  filters  are  called 
evanescent  mode  ( or  'ghost  mode'  as  Jaynes  called  it )  filters.  Microwave  bandpass 
filters  using  evanescent  modes  have  been  designed  successfully  [3,4,5].  The 
evanescent  mode  bandpass  filter  has  several  advantages  compared  to  the  conventional 

type  of  bandpass  filters  (waveguide  above  cutoff,  coaxial  line  etc.).  For  instance,  a 

sharper  transition  to  out-of-band  rejection  can  be  obtained  on  the  higher  frequency 
side.  Evanescent  mode  waveguide  filters  are  also  smaller  than  traditional  waveguide 
filters.  A  waveguide  operating  below  its  cutoff  frequency  is  basically  an  inductive 
element  [5,6].  Suitable  capacitive  elements  are  needed  to  construct  an  evanescent 
mode  waveguide  bandpass  filter.  Evanescent  mode  filters  using  conventional 
capacitive  elements  such  as  tuning  screws  [4,5]  are  costly  and  difficult  to  mass 
produce  because  of  their  complicated  structure.  Non-touching  E-plane  fins[7],  which 
are  easily  to  fabricated,  are  proposed  in  this  study  as  the  capacitive  elements.  The 
filter  structure,  shown  in  Fig.  II- 1,  consists  of  a  number  of  non-touching  E-plane  fins 
placed  in  a  rectangular  waveguide  below  cutoff.  The  fins  may  be  metal  only  or 
supported  by  a  dielectric  layer.  The  input  and  output  portions  of  the  filter  are  coupled 
to  the  external  circuits  via  double  step  waveguide  junctions.  The  larger  waveguides 
operate  above  the  cutoff  frequency.  Since  the  capacitive  elements  utilize  a  printed 
circuit  structure,  they  are  suitable  for  mass  production  at  low  cost. 

The  initial  approach  to  the  synthesis  of  an  evanescent  mode  filter  was  based 
on  image-parameter  theory  [3],  This  approach  cannot  be  used  for  exact  design 


technique  to  obtain  the  scattering  matrix  description  of  the  evanescent  mode  filter. 
The  insertion  loss  and  the  return  loss  of  the  filter  then  can  be  obtained  from  the  final 
scattering  matrix. 

Filters  designed  with  the  present  technique  have  been  tested  in  Ka  band 
(26.5-40.0  GHz).  Good  agreement  between  theory  and  measurement  is  observed. 


waveguide 
below  cutoff 


E-plane  fins 


double -step 
junction 


(a)  Cross  sectional  view;  (b)  Top  view; 


(c)  Side  view  of  the  filter 


II. 


ANALYSIS  AND  DESIGN  PROCEDURE 


The  analysis  and  design  of  an  evanescent  mode  waveguide  bandpass  filter 
with  non-touching  E-plane  fin  is  based  on  the  generalized  scattering  matrix  technique 
in  conjunction  with  the  spectral  domain  approach  and  the  mode  matching  method. 
The  study  begin  with  brief  description  of  the  generalized  scattering  matrix  .  Then  the 
scattering  matrix  representations  for  the  double-step  junction,  the  non-touching 
E-plane  fins,  and  the  evanescent  mode  waveguide  section  are  obtained  by  the  mode 
matching  technique,  spectral  domain  method,  and  waveguide  theory,  respectively. 
Finally,  these  scattering  matrices  are  combined  to  obtain  the  final  generalized 
scattering  matrix.  The  filter  response  is  calculated  from  the  final  scattering  matrix. 

A.  GENERALIZED  SCATTERING  MATRIX 

The  concept  of  a  generalized  scattering  matrix,  introduced  by  Mittra  and  Pace 
[10],  is  closely  related  to  the  scattering  matrix  of  circuit  theory  or  of  microwave 
network  theory.  It  differs  by  including  scattering  of  all  modes,  so  that  the  scattering 
matrix  will  in  general  be  of  infinite  order.  In  the  following  discussion,  the  term  " 
scattering  matrix  "  will  mean  "generalized  scattering  matrix  ." 

The  scattering  matrix  can  be  defined  for  a  junction  discontinuity  at  which  the 
fields  may  be  expanded  in  modes,  such  as  the  double-step  waveguide  junction  show 

in  Fig.II-3.  Consider  that  the  mnth  TE  mode  is  incident  upon  the  plane  z=0  from  the 

larger  waveguide  I  and  waves  are  reflected  back  into  waveguide  I  and  transmitted 
into  smaller  waveguide  II.  If  the  amplitude  of  the  mnh  TE  incident  in  I  is  normalized 


to  unity,  then  the  amplitude  of  the  pqth  TE  scattered  mode  in  waveguide  I  is  S  “ 
(pq.mn),  and  the  amplitude  of  the  pqth  transmitted  TM  mode  in  waveguide  II  isS“ 

(pq,mn)  .  The  terms  (pq.mn)  and  (pq.mn)  are  the  generalized  matrix 

elements  of  and  S^,  respectively.  The  other  scattering  matrixes  ,  such  as 
SEM  SEE  and  SMM 

,  etc.,  may  be  similarly  defined  [10,11].  In  conventional 
scattering  matrix  formulations,  the  propagating  modes  are  normalized  so  that  a  mode 
carries  unit  power.  Since  the  generalized  scattering  matrix  includes  evanescent 
modes,  such  a  normalization  is  inappropriate.  One  consequence  is  that  the  scattering 
matrix  is  not  symmetric. 

In  the  notation  used  above,  the  scattering  matrix  relates  the  excited  modes  to 
the  incident  ones  via : 


As 

=  [S] 

A*  ' 

_ 

SM 

si  ii 

A* 

Bs 

siii 

sii  ii 

B * 

The  superscript  s  indicates  scattered  fields,  and  i  expresses  the  incident  fields.  The 
general  element  of  S  is  (pq,mn),  where  x  and  y  =  E  or  M,  represent  a  TE  or  TM 
to  x  wave;  i  and  j  =  I  or  II  indicate  the  Region  I  or  Region  II;  m,  n,  p  and  q  are 
integers  corresponding  to  different  modes.  For  example,  Sjj^(ll.lO)  repressents 

the  ratio  of  the  amplitude  of  reflected  TMjj  in  Region  II  to  the  amplitude  of  an 
incident  TEjg  wave  in  Region  I.  In  more  detail  the  above  equation  runs  as  follows: 
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ASe(10)|  fSEE(10,10)  See(10,20)  ...  S^(10,01)  SeiJ|(10,02)... 

Ag(20)  See(20,10)  S^(20,20)  ...  SE^(20,01)  SE^(20,02)... 

A5Mm  S^(Ol.lO)  S^(01(20)  ..  SMM(01,01)  sMM(0i(02).. 

A^(02)  SWf(02,10)  S^(02,20)  ...S^(02,01)  S^(02,02). 

Bg(10)  SEE(10,10)  SEE(  10,20)  ...  Sf/fjdO.Ol)  SeM(10,02)... 

B|(20)  See(20,10)  See(20,20)  ...SeW(20,01)  SE1$20,02)... 


’a‘(io; 

Ag(20) 

A^(01) 

Ajy[(02) 

Bg(10) 

BE(20) 


bm(01)  SME(01,10)  S^(01,20)  ...SMM(0i,01)  sMM(0i>02)...  B>1) 

B^(02)  S^(02,10)  S^(02,20)  ...SMM(02,01)  S^(02,02)...  B^(02) 


(  II-2 ) 


Theoretically  the  generalized  matrix  is  of  infinite  dimensions  corresponding 
to  the  infinite  number  of  eigenmodes.  The  matrix  is  truncated  to  finite  size  for 
numerical  calculations. 

The  generalized  scattering  matrix  technique  is  useful  for  dealing  with 
structures  containing  several  discontinuities.  For  example,  the  structure  of 
evanescent  mode  waveguide  filter  with  non-touching  fins  consists  of  three  portions: 
the  double-step  waveguide  junctions  and  the  non-touching  E-plane  fin  portion,  and 
the  waveguide  below  cutoff.  Each  portion  of  the  structure  can  be  represented  by  its 
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corresponding  scattering  matrix.  These  matrices  are  then  combined  to  obtain  the 
scattering  matrix  of  the  filter  structure.  The  insertion  loss  and  the  return  loss  can  be 
obtained  from  the  related  elements  of  final  scattering  matrix. 


B.  SCATTERING  MATRIX  REPRESENTATION  OF 
A  DOUBLE-STEP  JUNCTION 

In  this  section  the  scattering  matrix  representation  of  a  double-step  junction 
is  derived  via  the  mode  matching  technique[12,13].  As  shown  in  Fig.  II-3  Region  I 
represents  the  larger  guide  and  the  smaller  guide  is  denoted  by  Region  II.  The  field 
is  derived  from  the  electric  vector  potential  F  and  the  magnetic  vector  potential  A 

[14]. 


E  =  -  VxF  +  Vx  Vx  A/j(oe0  (n-3) 

H  =  VxA  +  VxVxF/j  (op.0  (  H-4  ) 

In  this  study,  it  is  convenient  to  choose 

A  =  V  x  (  n-5  ) 

F  =  <j>  x  (  n-6  ) 


where  \|/  and  <j>  are  the  scalar  functions  which  represent  the  electric  (TE)  wave  and 
magnetic  (TM)  wave  transverse  to  the  x-direction,  respectively,  and  x  is  a  unit 
vector  in  the  x-direction.  These  potentials  are  extended  in  terms  of  their 
eigenfunctions  which  must  satisfy  the  boundary  conditions.  The  potentials  in  each 


region  can  be  expressed  as  follows: 


M  N-l 

S  l  exP(  'kizhmnz  >  +  bimn  exP(  ^hmn2  (  n‘7  ) 

m=l  n=0 

M-l  N 

—  ^  ^  f  cimn  exP(  "kizemnz  )  +  ^imn  exP(  kizemnz  ^Qimn^x,y^  (  ) 

m=0  n=l 

( The  notation  used  in  this  part  is  described  below: 

X  i  c  j  d  1 

1  =  m,  n  or  mn  indicates  mth,  nth  or  mnth  mode 
d  =  h  or  e,  h  :  TE-to-x  field  ;  e  :  TM  -  to  -x  field 
j  =  1,  2,  or  3  indicates  the  variable  in  the  region  1,  2,  or  3  of 
the  smaller  waveguide. 

c  =  x,  y,  or  z  indicated  the  x- ,  y-  or  z-direction,  respectively 
i  =  I  or  II,  I :  larger  waveguide;  II :  smaller  waveguide 
indicate  the  variable ) 

where 

p.mn(x>y)  =  Rim(x)Sin(y)  (  n-9  ) 

Qimn(x-y)  =  TimWVin(y)  (  D-10 ) 

are  the  eigenfunctions  of  the  TE-to  -x  and  TM-to-x  field,  respectively.  In  region  I, 
these  eigenfunctions  are  represented  by  : 


RImW  =  Nrlm  sin  ( m7lx  1  ai> 
Sin(y)  =  NSIn  cos(  ntty  /  bj) 


(0-11) 
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(0-12) 


Nnm-V^TiJTI,  NSIn=V(8n0+ 1)  /  b  <n-13) 

m=l,2, 3,  —  M  n  =  0, 1, 2,  •••  N-l 

TlmW  =  Nqjjj  cos  ( mroc  /&{)  (  n- 14 ) 

vm(y)  =  NvInsin(n7Iy''bl'  (  O'15 ) 

NtIm  =  ''<«mo+1>/al  Nyln  =  V  (8„0- 1)  /  aj  (n-16) 

m  =0, 1, 2,  -  M-l  n  =  1, 2,  3,  -  N 

where  8mo  and  8no  are  the  Kronecker  delta  functions.  In  Region  II,  the 
eigenfunctions  are  given  by  : 

snq(y)  =  Nsnq  cosK  W  /  bi)(y-yi)3  ( 0-17 ) 

Nsnq=''<V1)/b>  (n'18) 


Xlllp  sinknxlhp(x-Xl) 

NrH2p  sin  kIIx2hp(  x_xrhl) 

Xj  <x  <Xj  +  hj 

RnPW=< 

+  Nrnmpcos  krn2hp(  x’xrhi) 

X1  +  bl  <  x  <  xl  +  bl+  b2 

\  Nrii3p  sinkIIx3hp(  x-x-a) 

Xj  +  h|  +  h2<  x  <  Xj  +  a 

(  0-19) 

p  =1,  1,  2, 

P  q  =  0,  1,2,  *•*  Q-l 

vnq(y)  =  Nvnqsin  K  <1*  /  bi)(y-yi)] 

( 11-20 ) 

^vOa  =  ^  ^no‘  /  al 

(  0-21 ) 

Tllp(x)  - 


Ntnip  °°s  knxiep(  x'xi' 
^tll2p  008  knx2ep(  x'xrhl) 

+  Ntnp  sln  kiix2ep(  x'xrhi) 
>NtII3psinkIIx3ep<  x'xra) 


Xj  <  x  <  Xj  +  hj 


Xj  +  hj  <  x  <  Xj  +  hj+  h2 


Xj  +  hj  +  h2<  x  <  Xj  +  a 


(  11-22 ) 


p  =0,  1,  2,  ••  P-1  q  =  1,  2,  3,  -  Q 


where  NrIIj  and  N^j j  =  1, 2,  and  3,  are  normalized  coefficients  of  eigenfunctions 

Rnm  (x)  and  Tjj^x)  of  the  smaller  waveguide  in  region  j,  respectively.  They  are 
given  in  Appendix  C.  The  mth  eigenmodes  of  TE-to-x  and  TM-to-x  fields  in  the 
partially  filled  waveguide  in  region  j  are  represented  by  knxjhm  311(1  kHxjem  (They 


correspond  to  Yim  and  Yim  in  Part  one  of  this  study  and  can  be  found  as  described 

in  Part  one).  Here  and  kjzemn  are  the  propagation  constants  of  the  mn  th  TE 

and  TM-to-x  fields  in  the  z-direction,  and  must  satisfy  the  following  dispersion 
equations 


k2lzdmn  +  k02  =  ( mtc  /  a^2  +  ( n nl  bj)2 
k2llzdmn  +  kj2  =  k2llxjdm  +  ( nn  /  b)2 


(  H-23 ) 


(  H-24) 


where  d  =  h  or  e,  j  =  1,  2,  and  3,  kQ  is  the  wave  number  in  free  space,  kj  =  Vijk0, 

and  £■  is  the  relative  dielectric  constant  of  the  smaller  guide  in  region  j.  The 

coefficients  a^,  bimn,  cimn,  and  dimn  in  Eqs.  ( II-7  )  and  ( II-8  )  correspond  to 

incident  and  reflected  waves  and  are  related  to  each  other  by  the  scattering  matrix. 
The  scattering  matrix  can  be  determined  by  matching  the  tangential  fields  at  the  step 
discontinuity  at  z  =  0.  ( Fig,  El-4  ) 

^Ix,y  I  at  z=o  =  Bnx,y  I  at  z=o  3163  ^2 

^lx,y  I  at  z=o  =  ^llx,y  lat  z=o  *n  area  ^2 

^lx,y  I  at  z=o  =  0  area  -  A2  ( D-25  ) 


Applying  the  orthogonality  relationship  between  potential  functions  leads  to  the 


matrix  equation : 


Maj  0  -V5n  0  ' 

^imn 

0  Vd  0  -Mdll 

^•imn 

Mbl  Md  -VaII  Vdn 

Bnpq 

val  ‘vdl  ^bn  Mcn  t 

DIIpq 

Mi 


■Mai  0  Vbn  0 

Blmn 

0  -VCI  0  MdII 

®Imn 

-Mbi  Mcl  Van  Vdn 

^Ilpq 

.  val  vdl  Mbll  ‘Men 

v - v - /  ' 

cnpq 

t 

M2 


where  Maj)  Mcj...ect.  express  matrices.  For  instance,  is  a  matrix  of  dimensions 
PQ  x  MN.  Vjju  Vaj ...  are  diagonal  matrix  of  dimensions  PQ  x  PQ,  MN  x  MN  ..., 


'Vi 


|  tUt'i  »*l  i"l.k>l  ■*>  t’i  Vt.t'L 
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■^mn>  Dpq-  .cct  are  vectors  of  dimensions  MN  or  PQ,  and  are  given  in  Appendix  D. 

The  scattering  matrix  of  the  double-step  discontinuity  is  then  given  by 

f  N. 

Sdll  Sdl2 

Sd=  = 

s  S(J22  (  D-26 ) 

To  check  the  validity  of  the  calculations,  we  calculated  the  scattering 
parameters  of  a  double-step  waveguide  discontinuity  from  Ku  to  K  band  ( 15.8  x  7.9 

mrn^  to  10.7  x  4.32  mm^ )  and  compared  results  with  Arndt  and  Wriedt  ’s[15].  The 
results  are  shown  in  Fig.  11-5.  They  are  in  good  agreement .  Also,  the  magnitudes 
of  the  scattering  parameters  of  resonant  irises  with  finite  thickness  t  ( Fig.  H-6 )  are 
obtained  by  combining  two  step  junctions  with  a  evanescent  mode  waveguide  of 
length  of  t.  The  numerical  results  are  compared  with  available  data  [12]  and  are 
shown  in  Fig.II-7  and  Fig.  II-8.  Once  again  close  agreement  is  observed. 


SCATTERING  MATRIX  REPRESENTATION  FOR  THE  PORTION 
OF  NON-TOUCHING  FINS 


The  scattering  matrix  for  a  non-touching  fin  can  be  obtained  by  combining 
the  spectral  domain  method  with  residue  calculus.  ( The  details  are  described  in  the 
first  part  of  the  dissertation.)  Since  the  filter  may  contain  more  then  one  fin  element, 
a  general  case  must  be  considered  here.  Figure  II-9  shows  that  a  structure  consists  of 

n  fins  ,  which  are  indicated  by  Fj,  F2, ...  and  Fn.  These  fins  are  separated  from  one 
to  another  by  distances  Tj,T2,—  and  Tn_j.  To  obtain  a  scattering  matrix 
representation  Sf  for  this  portion,  we  consider  that  a  substructure  which  consists  of 
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Present  method 
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Frequency  (GHz) 


Reflection  coefficient  I  Sjjl  as  a  function  of  frequency,  if  a  TEjq 

wave  is  incident  from  Ku-band  waveguide  ( 15.8  x  7.9  mm2  )  to 
K-band  guide  ( 10.7  x  4.32  mm2 ) 


Frequency  (GHZ) 


Fig.  n-7  Magnitude  of  reflection  coefficient  of  resonant  iris  with  finite  thickness  t 


aj  =  15.8  mm,  bj  =  7.9  mm,  a  =  1 1.2  mm,  b  =  5.6  mm,  t  =  2  mm 
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Fig.  II-8  Magnitude  of  transmission  coefficient  of  resonant  iris  with  finite 
thickness  t.  aj  =  15.8  mm,  bj  =  7.9  mm,  a  =  1 1.2  mm, 
b  =  5.6  mm,  t  =  2  mm 
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fins  Fj  and  F2,  with  spacing  Tj.  The  substructure  is  shown  in  Fig.  II- 10.  The 

scattering  matrixes  of  fins  Fj  and  F2,  Sfj  and  Sf2,  are  obtained  by  the  method  given 

in  Part  one.  With  the  scattering  matrix  notation  described  above,  each  of  these 
matrices  contains  four  submatrices  and  can  be  written  as 


sfl 


sf2 


Sfm 

Sfll2 

Sfl21 

Sfl22 

Sf2ll 

Sf2l2 

Sf221 

Sf222 

( n-27 ) 


( n-28 ) 


With  the  knowledge  of  the  scattering  parameters  for  a  single  fin,  the 
generalized  scattering  matrix  technique  is  applied  to  obtained  a  scattering  matrix  of 
this  substructure.  The  concept  is  that  of  multible-reflection  phenomenon.  If  a  wave 

from  Region  A  is  incident  upon  Fj,  fields  will  be  partly  reflected  back  into  Region  A 

and  partly  transmitted  into  Region  B.  After  traveling  a  distance  Tj,  a  part  of  the  wave 
transmitted  into  Region  B  is  reflected  back  and  the  remaining  part  is  transmitted  into 
Region  C  via  F2.  This  process  continues  until  the  reflected  wave  dies  out.  This 

multiple-reflection  phenomenon  between  the  Fj  and  F2  implies  a  matrix  combination 
process  that  leads  to  the  scattering  matrix  for  the  substructure. 

In  Region  B  of  Fig.  II- 10,  the  smaller  waveguide  section  with  length  Tj, 


operated  below  its  cutoff  frequency,  is  represented  by  the  matrix  ST1.  The  elements 


of  those  matrices  are  easily  obtained  from  the  transmission  line  equivalent  circuits  of 
the  waveguide,  in  which  the  characteristic  impedances  are  imaginary.  Since  no 
propagation  modes  exist  below  cutoff  frequency,  this  waveguide  section  in  the 
equivalent  circuit  acts  as  a  lumped  reactance.  The  wave  travels  a  distance  T1  in  this 

guide  so  each  mode  is  multiplied  by  a  exponential  decay  factor  exp(-KjIwj  Tj), 


where  d  =  h  or  e  represent  TE  or  TM  waves,  respectively.  In  more  detail ,  ST1  can 


be  written  as 

ST1  = 


where  the  submatrix  Stjis  given  by 


( 11-29 ) 


exp(-  K{IzhioTl)  0  .  .  .  0 

0  exP('  Kllzh20Tl^ 

exp(-  Kjj^qjTj) 

0 

exP(-  KHzepQTl) 


( H-30) 


The  combination  of  Stlwith  Sfjand  Sf2  results  in  the  scattering  matrix  Sj 


in  Fig.  II-9  that  represents  the  cascaded  structure.  The  algebraic  process  is  shown  in 
Appendix  E.  The  elements  of  matrix  S,  are  given  by 


D.  SCATTERING  MATRIX  REPRESENTATION  OF  CASCADED 

SECTIONS 


The  side  view  of  the  filter  can  be  now  represented  by  Fig.II-11.  The 
double-step  junctions  are  described  by  scattering  matrices  S^j  and  S^- The 
capacitive  element  of  the  filter  is  the  portion  of  the  non-touching  E-plane  fins,  which 
is  represented  by  scattering  matrices  Sf.  The  matrix  ST  represents  the  evanescent 

mode  guide  section,  and  it  can  be  found  in  the  same  way  as  ST].  Since  S^j, 


4<»  .»*  I  ,  i  .  t  J  «  J.t  M  4.*  M  l.f  f.i  M  * J  V Vl I**  V  i/fc 


SAA  =  Sflll+  Sf  1 12  Stl  Sf2l  1  Stl  Sfl21 

SAC  =  Sf  1 12  Stl  U2  Sf2l2 

SCA  =  Sf221  Stl  U1  Sfl21 

SCC  =  Sf222  +  Sf221  Stl  U1  Sfl22  Stl  Sf2l2 

^1  =  (  I  -  Sfl22^tl  Sf2ll  Stl)_1 
U2  =  (  I*Sf2nStl  Sfl22Stl)‘1 

where  I  is  the  identity  matrix.  In  the  next  step,  we  consider  the  substructure  that 
consists  of  fins  Fj,  F2,  and  F3  with  spacing  Tj  and  T2.  The  combination  of  Flt 

F2,  and  Tj  is  now  expressed  by  matrix  Sj.  Using  ST2  and  Sf3  to  represent 
scattering  matrices  of  T2  and  F3  respectively,  we  obtain  the  matrix  S2  of  this 
substructure  in  the  same  way  as  that  of  deriving  matrix  Sj.  The  same  procedure  is 
repeated  until  the  matrix  Sf  is  obtained  ( Fig.  II- 9). 


P 


Sf,  and  ST  are  all  known,  the  overall  matrix  S,  which  characterizes  the  filter 

structure,  can  be  determined  from  these  scattering  matrices  by  the  generalized 
scattering  matrix  technique.  The  elements  of  the  matrix  S  are  given  by  : 


rO) 

/  N 

S 1 1  Si2 

/  \ 
10) 

r(2) 

S21  S22 

1(2) 

k 

/ 

where  the  submatrices  Sjj,  S12,  S21,  and  S22  are  repraesented  by: 

S11  =Sdlll  +Sdll2SD"lsu 
^12  =  ^dll2+  SD_1SV 
®2l  =  ^dll2  + 

S22  =  Sdlll  +Sdll2SD'lsQ 

SD  =  ( I  -  St  Sfl  j  St  Sdl22  )  ( I  -  St  Sft2  StSd2ll  )  - 

StSfl2  StSd2ll  StSf21  StSdl22 

SU  =  StSfU  st  Sdl21  ( I  -  StSf22  St  Sd2ll  )  + 

St  Sfl2 St  Sd2ll  St  Sf2i  St  Sdl2l 
SV  =  St  Sfl2  StSd2l2  ( I  -  St  Sf22  St  Sd2ll  )  + 

StSfi2StSd2ii  St  Sft2  St  Sd2l2 
SP  =  ( I  -  StSfn  StSdj22  )  St Sf2!  stsdi21  + 

®t  Sfi i  StSdi2i  stsf21  stsdl22 

SQ  =  ( I  -  st  Sfll  St  sdl22 )  St  Sf22  St  Sd2l2  + 

stsfi2  stsd2i2  stsf2i  stsdl22 
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DESIGN  EXAMPLES  AND  DISCUSSIONS 


The  behavior  of  a  waveguide  below  cutoff  is  basically  that  of  an  inductive 
element.  Capacitive  elements  are  needed  for  a  bandpass  filter.  The  non-touching 
E-plane  fin  may  be  represented  by  its  equivalent  r-network  (  Fig.  1-2  (d)).  The 
element  values  in  the  equivalent  network  can  be  found  from  Sf.  In  principle,  the 

width  and  height  of  the  fin  have  to  be  chosen  such  that  X2  ( in  Fig.  1-2  (d))  is 
negative,  equivalent  to  a  shunt  capacitor. 

Several  filters  have  been  designed  by  the  procedure  described  above.  The 
design  procedure  can  be  used  in  a  wide  frequency  range.  Howevre  in  this  study, 
only  bandpass  filters  operating  in  Ka  band  have  been  considered.  In  the  following 

discussion,  the  larger  waveguide  is  WR-28  (7.1 1  mm  x  3.56  mm)  while  the  smaller 

waveguide  is  WR-15  (  3.76  mm  x  1.88  mm).  The  center  frequency  is  the  most 
important  quantity  for  design  of  a  bandpass  filter.  Attention  is  first  directed  to  how 
the  height  and  width  of  the  single  fin,  and  the  distance  between  the  edge  of  the  fin 
and  the  step  junction  affect  the  center  frequency  of  the  evanescent  filter. 

Figure.  H-12  shows  the  relationship  between  the  height  of  the  fin  and  the 
center  frequency  for  a  filter  with  the  E-plane  fin  supported  by  a  Duroid  substrate 

(dielectric  constant  e  =  2.2,  thickness  h2  =  0.127  mm).  The  center  frequency 
decreases  as  the  height  d  of  the  fin  increases,  because  as  d  increases  there  is  more 
electric  energy  stored  in  the  gap  between  the  fin  and  the  wall  of  the  guide.  This 
corresponds  to  a  larger  shunt  capacitor  in  its  equivalent  network  while  the  series 
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inductance  Xj  exhibits  little  change. 

Figures  11-13  (a)  and  (b)  show  the  relationship  between  the  width  of  the  fin 
and  the  center  frequency  for  a  fin  with  a  fixed  height.  The  wider  the  fin,  the  lower 
the  center  frequency  is.  The  wider  fin  leads  to  a  larger  capacitance  and  inductance  in 
its  equivalent  circuit  and  hence  a  lower  resonant  frequency.  (  See  Fig.  1-12) 

The  insertion  loss  versus  the  frequency  for  different  distance  T  of  the  filter 
with  fixed  width  and  height  of  the  fin  is  shown  in  Figs.  11-14  and  11-15.  In  Fig. 
11-14  the  fin  is  supported  by  a  Duroid  substrate  while  in  Fig.  11-15  there  is  no 
dielectric  substrate.  When  the  fin  is  separated  away  from  the  step  junction,  the 
center  frequency  becomes  higher  and  the  curve  becomes  steeper.  This  phenomenon 
can  be  explained  as  following.  To  simplify  our  discussing,  we  assume  that  only  the 

TE10  evanescent  mode  exists  in  the  smaller  guide  containing  a  narrow  fin  without 
substrate.  The  waveguide  can  be  described  by  a  simple  transmission  line  equivalent 
circuit  with  imaginary  characteristic  impedance  ZD  =  jx0  (Fig.  11-16  (a)),  or  in  terms 

of  its  7t-equivalent  network  as  shown  in  Fig.  11-16  (b).  The  fin  is  represented  by 
the  equivalent  T-network  (Fig.  11-16  (c)).  If  the  value  of  X]  is  small  compared  to  X0, 

then  the  equivalent  circuit  of  the  filter  may  be  expressed  by  Fig.  11-16  (d),  where  X0 
is  given  by 

X0=  I20brc/M(XAC)5-1 )  (11-32) 

The  propagation  canstant  of  the  TE10  mode,  y,  is  given  by 

y=  2tcV  (XAc)2-1  A  (n-33) 
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in  which  X  is  the  free-space  wavelength  and  Xc  is  the  cutoff  wavelength. 

The  resonant  frequency  of  the  filter  in  Fig.  11-16  (d)  is  mainly  determined 

by  X2  and  the  parallel  inductance  Lp=  jX0coth  (0.5yT).  When  T  becomes  larger,  Lp 
becomes  smaller  while  X2  is  unchanged.  Hence  the  resonant  frequency  becomes 

higher.  If  T  is  large  enough,  Lp  approaches  to  jXQ  since  the  term  of  coth(0.5yr) 
tends  to  1.  Hence  for  large  T,  the  center  frequency  of  the  filter  approaches  to  a 
certain  value.  Also  the  curve  becomes  sharper  when  T  increases,  since  the  coupling 
between  the  fin  and  the  step  is  weaker. 

With  the  data  presented  above,  it  is  now  possiple  to  design  a  filter  with  one 
fin  element.  Fig.  11-17  shows  the  results  for  a  filter  designed  in  Ka  band  using  one 
fin  element  without  a  dielectric  layer.  The  solid  curve  represents  the  results  obtained 
by  this  analysis  and  the  dashed  line  indicates  the  measured  data.  They  are  in  good 
agreement.  Figure  11-18  shows  the  results  of  the  filter  with  the  E-plane  fin 
supported  by  a  Duroid  layer.  Once  again  agreement  between  theoretical  prediction 
and  experimental  data  is  quite  good.  The  small  deviation  of  the  insertion  loss 
between  theory  and  experiment  at  the  center  frequency  comes  from  the  metal  and 
dielectric  loss. 

Figure  11-19  shows  the  calculated  response  of  a  filter  that  consists  of  two 
equal  E-plane  fin  elements  on  a  Duriod  substrate,  The  height  of  the  fin,  the  distance 
from  the  edge  of  the  fin  to  the  double- step  junction,  and  the  spacing  of  the  two  fins 
are  the  design  parameters.  The  3  dB  bandwidth  is  about  1.6  GHz.  It  is  noted  that 
the  steepness  of  the  out-of-band  insertion  loss  curve  on  the  higher  frequency  side  in 


Fig.  11-19  is  almost  equal  to  the  one  or  *he  lower  frequency  side.  This  happens 
because  the  waveguide  below  its  cutoff  frequency  acts  like  a  lumped  reactance. 

In  general,  wider  bandwidth  and  better  transmission  performance  can  be 
achieved  by  increasing  the  number  of  fins.  For  a  multi-fin  structure  an  optimization 
procedure  similar  to  [161  may  be  used  to  optimize  the  performance  of  the  filter. 
Alternatively,  a  filter  synthesis  can  be  used  to  find  the  required  equivalent  circuit 
parameters  in  the  filter.  The  necessery  fin  dimensions  and  fin  spacing  can  be  found 
from  a  look-up  table.  The  center  frequency  can  be  controlled  by  the  dimensions  of 
the  fins  and  the  distance  between  edge  of  the  fin  and  the  double-step  junction.  The 
shorter  the  height  and  the  narrower  the  width,  the  higher  the  center  frequency.  Also 
a  longer  distance  from  the  edge  of  the  fin  to  the  double-step  junction  leads  to  a 
higher  central  frequency.  The  height  of  the  fin  is  the  most  sensitive  parameter  to  the 
center  frequency.  For  instance,  in  Figs.  II- 13  (a)  and  (b),  if  the  width  of  the  fin 
chages  by  0.1  mm,  the  certer  frequency  shifts  about  0.5  GHz,  while  the  variation 
of  0.05  mm  on  the  height  corresponts  to  about  1  GHz  frequency  changing  (See 
Figs.  IM2  (a)  and  (b)). 

Since  the  non-touching  fins  can  be  produced  by  photolithographic 
techniques,  fine  tuning  is  normally  not  required. 


Frequency  (GHz) 


Fig.  11-12  (a)  Return  loss  of  the  filter  versus  frequency  for  different  height  d  of 

the  E-plane  fin.  hi  =  1.753  mm,  h2  =0.127  mm,  ^  =  2.2  w  =  0.3  mm 

T  =  2.0  mm,  A  :  d  =  1.4  nun.  B:  d  =1.45  mm,  C  :  d  =  1.50  mm, 

D  :  d  =  1.55  mm,  E  :  d  =  1.60  mm 


Frequency  (GHz) 


Fig.  II- 12  (b)  Insertion  loss  of  the  filter  versus  frequency  for  different  height  d  of 

the  E-plane  fin.  hi  =  L753  mm,  h2  =0. 127  mm,  =  2.2  w  =  0.3  mm, 

T  =  2.0 mm,  A:d=1.4mm,  B:  d  =1.45  mm 
C  :  d  =  1 .50  mm,  D  :  d  =  1 .55  mm,  E  :  d  =  1 .60  mm 


Frequency  (GHZ) 


Fig.  11-13  (a)  Return  loss  versus  frequency  for  different  width  of  E-plane  fin 
d  =  1.483mm,  T  =  2.0mm,  A  :  w  =  0.1mm,  B:  w  =0.2mm 
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Fig.  II- 1 3  (b)  Insertion  loss  of  the  filter  versus  frequency  for  different  width  w  of 

E-plane  fin.  hi  =  1.753  mm,  h2  =0.127  mm.  e2  =  2.2,  d  =  1.483  mm, 

T  =  2.0  mm,  A:w  =  0.1mm,  B:  w  =0.2  mm 
C  :  w  =  0.3  mm,  D  :  w  =  0.4  mm,  E  :  w  =  0.5  mm 
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Fig.  H-15  Insertion  loss  of  the  filter  versus  frequency  for  different  distance  T 
between  the  step  junction  and  the  E-plane  fin  (metal  only), 
d  =  1.50  mm,  w  =  0.28  mm,  A  :  T  =  1.5  mm,  B:  T  =  2.0  mm 
C  :  T  =  2.5  mm,  D  :  T  =  3.0  mm,  E  :  T  =  4.0  mm 


(a)  jx0coth(YT/2) 

(b) 


jx0coth(YT/2) 

(d) 


Fig.  11-16 

(a)  Transmission  line  equivalent  circiut  for  evanescent  mode  waveguide 

(b)  Tt-network  equivalent  circiut  for  evanescent  mode  waveguide 


(c)  T-network  equivalent  circiut  for  metal  fin 

(d) A  simplified  euqivalent  circuit  of  the  evanescent  mode  waveguide  filter 


Frequency  (GHz) 


Fig.  II- 17  Computed  and  measured  insertion  loss  for  an  evanescent  mode 

waveguide  filter  with  an  E-plane  fin.  (metal  only,  =  1.753  nm 

h2  =  0.157  mm,  d  =  1.511  mm,  w  =  0.787  mm,  T  =  2.506  mm) 
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Fig.  H-18  Computed  and  measured  insertion  loss  for  an  evanescent  mode  waveguide 
filter  with  an  E-plane  fin  on  a  Duriod  substrate.  (hj=  1.753  mm,  £2  =2.2 
h2  =  0.157  mm,  d  =  1.483  mm,  w  =  0.292  mm,  T  =  1.502  mm) 
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Fig.  11-19  (a)  Insertion  loss  for  an  evanescent  mode  waveguide  filter  with  two 
E-plane  fins  on  a  Duriod  substrate.  (hj=  1.753  mm,  h2  =  0.127  mm. 

d  =  1.50  mm,  w  =  0.292  mm,  T  =  2.0  mm,  Tl=  4.6248  mm,  e-,  =  2.: 
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Fig.  H-19  (b)  Return  loss  for  an  evanescent  mode  waveguide  filter  with  two 

E-plane  fins  on  a  Duriod  substrate.  ( hj=  1.753  mm,  h2  =  0.127  mm, 

d  =  1.50  mm,  w  =  0.292  mm,  T  =  2.0  mm,  Tl=  4.6248  mm,  £2  =  2.2) 
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IV  CONCLUSIONS 


The  analysis  and  design  of  evanescent  mode  waveguide  filters  with 
non-touching  E-plane  fins  are  presented.  In  the  analysis,  the  fundamental  mode  as 
well  as  the  higher  order  modes  effects  have  been  taken  into  account.  This  is 
important  in  an  accurate  filter  design.  The  filter  designed  in  Ka  band  has  been  tested 
and  good  agreement  between  the  measurement  and  the  theory  is  observed.  It  is 
believed  that  this  kind  of  filter  will  be  useful  in  the  millimeter-wave  frequency 
region. 
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APPENDIX  A.  CLOSED  FORM  EXPRESSION  OF  GREEN’S  FUNCTION 
IN  THE  SPECTRAL  DOMAIN 


According  to  [9],  the  TM-to-x  (LSM)  and  TE-to-x  (LSE)  equivalent 
transmission  lines  for  the  E-plane  strip  described  in  Fig.  2  can  be  drawn  in  Fig.  A-l. 
Here, 


Y t-jal+P'-i'k2 

y, 


"TMi 


2n, - » "1,2,3,  n  «  0,1,2, 3  •  •  • 

Y/ 


(  A-l  ) 
(  A- 2  ) 

(  A-3  ) 


The  driving  point  input  impedance  Ze  for  the  TM  equivalent  circuit  is  given  by 
1 

Z'« - 

yi  +  y{y 

where 

yiy  “  y-TMy^^y^i 

t  y~TMi  +  yiCothy;^; 

yi  y™2y{  +  yrM2colhy2h2 

Zi'-ZTMiCOthYiAi. 

Notice  that  ye  and  y *  are  input  admittances  looking  left  at  x  =  hj  and  x  =  hj  +  h2) 
respectively,  while  ye2ris  the  admittance  looking  right  at  x  =  hj  +  h2  Similar 


(  A- 6  ) 
(  A- 7  ) 


equations  can  be  written  for  the  TE  equivalent  circuit: 


Z*.  1 

A*  L  L 

(  A- 8  ) 

yi  +  yiy 

yir  “  J'TE3colhYjA, 

(  A-9  ) 

-  y-reiCO^Y!*! 

(  A- 10  ) 

.  -Kxe2  +  ^i*cothy2*2 

yi-y-n 1  yf  +  yJZ1coihy7h1 

(  A- 1 1  ) 

Finally,  the  Green, s  functions  in  the  spectral  domain  can  be  represented  by  Z*1  and 

Zc  as  follows  [8]: 

Gyy  ■  Z'N}  +  ZkNy 

(  A- 12  ) 

Gy,-G,y-{Z'-Zh)N'Ny 

(  A- 13  ) 

G,,  “  ZkN?  +  Z'Ny 

(  A- 14  ) 

where 

(  A- 15  ) 

Note  that  the  denominators  of  Z*1  and  Ze  are  the  eigenvalue  equationn  of  LSE  and 

LSM  modes. 

Y,.?TMI  YJ.ZTM3  Y|  ?Ul 

Y*,7UZ  Ti.ZttS 

X=nt  xsllfft;  X=a  x-0 


x  =  h,  x  =  fii  *  ri 2  x  =  o 


Fig.  A-l  The  TM  and  TE  equivalent  transmission  lines  for  the  E-plane  strip 


APPENDIX  B.  EXPRESSION  OF  NORMALIZED  REACTANCE  Xj  AND  X2 


The  equivalent  circuit  element  values  of  the  non-touching  E-plane  fin  are 
derived  as  follows.  In  Fig.  B-l,  the  normalized  input  inpedance  Zin  can 

be  expressed  as 


Zin=jX1+Z" 

jX2(l+jX1) 

=  jXi+  - 

1  +  j  (Xj  +  X2  ) 

jXj-  X1(X1-X2)+  jXjd+jXj) 

1  +  j  (Xj  +  x2 ) 

1 

^in  +i  ®in 

where 

1  +  j  (Xj  +  X2 ) 

^in  +  J  ®in  = 

j  (Xj  -f  X2  )  -  Xj  (Xj  +  2X2) 


(B-  1) 


( B-2 ) 


Equation  (B-2)  can  be  repressed  as 

1  +  j  (Xj  +X2 )  =  j  Gjj,  (Xj  +  X2)  -  Bin  (Xj  +  X2  ) 

-  Gin  X,  (Xj  +  2X2 )  -  jBinXj  (X,  +  X2  )  (  B-  3  ) 

The  real  part  and  imaginary  part  of  the  left  hand  side  of  the  Eq.  (B-3) 
should  be  equal  to  those  of  the  right  hand  side,  a  set  of  coupled  equations  is  obtaind: 


,\  , 


,v 


c 

/ 


( B-4 ) 


Gin(X1  +  X2)  -  BinX1(X1  +  X2)  =  (X1  +  X2) 

Bin(Xl  +  X2)  +  GlnX,(X1  +  2X2)  =  -l  (  B-5  ) 


Xj  and  X2  are  found  by  solving  the  Eqs.  (B-4)  and  (B-5).  The  result  is  given  by: 


( B-6 ) 

( B-7 ) 


Jx,  jX, 


Zin 

Z,=  l 

)*! 

4 

Zin 

||  1*2 

Z1  = 

1+jXj 

JX, 

jX2(l+jX1) 

^  1 

Zin 

• 

Z|"  =  ' 

1  +j(X1+X2) 

Fig.  B-l  Equivalent  circuit  to  find  the  normalized  input  impedance. 
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APPENDIX  C.  NORMALIZED  COEFFICIENTS  OF  THE  EIGENFUNCTIONS 


IN  REGION  II 


In  reference  to  Fig.  II-4  and  Eqs.  II-3  and  II-4,  the  tangential  fields  in 


Region  II  are  given  by 


Eny  =  -  3<(>2  /  dz  +  (1  /  jcoe0£j)  3fy2 1  dxty 


-  -  Z  X  kuzhmn^Ilmn  +  ®IImn)  EIImn 
m=l  n=0 


(1/jco  £j  eQ)X  X  (nn/b)cos[n7t/b(y-y  j)]  (-CIImn  +  DIImn)  3Tm/  dx 
m=0  n=l 

(C-l) 


Eiiz  =  3<t>n  /  +  C1  /  jcoe0ej)  d^H  /  3xdz 


=  -  X  X  (nrc/b)sin[n7t/b(y-y j)]  (-AIImn  +  BIImn)  RIImn 


m=l n=0 


(l/jcoe0)  X  X  knzemn  (-CIImn+  DIImn)  n  sin[n7c/b(y-yi)]  aTm/3x  /  ej 


m=0  n=l 


( C-2 ) 


Hjjy  =  a^n  /  +  0  /  jOHi0)  >  ^x5y 


-  X  X  ^nzemn^^-'llmn  +  ^Ilrnn)  Qllmn 

m=0  n=l 


(1  /  jcojj.0)  X  X  (n7t/b)sin[n7t/b(y-y,)]  (AIImn+  BIlmn)aRm/  dx 


m=l  n=0 


^  .,na  a  v  v  *.•  vnj% 'S  ’A-  v  * »  * .  *.  A  » >* 


.%  .%  .%  > 


( C-3  ) 


Hnz  =  +  0  /  jo>p.0)  92<j>n  /  dxdz 

oo  oo 

=  - 1 1  (nJc/b)cos[n7t/b(y-y!)]  (CIImn  +  DIImn)  TIImn  + 
m=On=l 


(1  /j©*i  )I  (-AIImn  +  BIImn)  sin[n7t/b(y-yi)]  9Rm/3x 


where 


m=l  n=0 


(  C-4 ) 


^Ilmn  ~  a  IlmnexP  (~kIIzhmnz) »  ®nmn  ~  b  lImnexP  (kllzhmnz) 


^Ilmn  ~  cIlmnexP  (‘kIIzemnz) »  ^Ilmn  “  ^  IlmnexP  (kIIzemnz) 

(  C-5 ) 

The  following  equations  are  obtained  by  matching  the  tangential  fields  at  the 
interface  x  =  xl  +  hj  and  x  =  Xj  +  +  h2. 


Nrli2m  =  Nrllimsin  (kllxihm  hP 

kIIx2em  ^tl!2m  =  'kllxlem  ^tnims^n  (kllxiem  *h) 


(C-6) 


( C-7  ) 


e2Ntll2m  =  eiNtnimcos  (kllxlem  (  C‘^  ) 

kIIx2hm  ^rn2m  =  kIIxlhm  ^rlllmcos  ^kIIxlhm  (  C-9  ) 

'NrII3msm  (kIIx3hm  h3^ 

=  ^ffl2ms*n  (kIIx2hm  ^2)  +  ^rII2mcos  (kIIx2hm  ^2)  (  C-l( 

kIIx3em  NtII3msin  (kIlx3em  ^3 

=  ‘kIIx2em  P^tll2ms*n  (kllx2em  ^2)  "  ^tll2mcos  ^Ilx2em  (C-l] 


V.V_V.V.VvV.N' -.v-'.v '  -  V-  V* V-V- 


( C-12 ) 


'e3NtII3mC0S(kIIx3em  h3^ 


“  ^tII2mcos  (kllx2em  ^  +  ^tll2ms^n(kllx2em  ^2^ 


kIIx3hm  Nrn3mcos  (kIIx3hm  h3) 

=  kIIx2hm  t^rII2mC0S  (kIIx2hm  ^2)  ‘^rll2ms'n  (krix2em  ^2^  (  C- 13  ) 


The  coefficients  Nrnj  and  NUIj  are  normalized  such  that 


Xj  +  a  Xj  +  a 

JRnm(x)dx=!  and  J  TIIm(x)dx  =  l 


(  C-14) 


Noimalization  leads  to  the  results 


^rlllm  "  ^.5  (  *  +  ^  ^  0-5hj  -  0.25sin  (2  knxihm  ^iVkiixlhm  J  + 

l  kllx3hm  cos(  kllx3hmh3)i2 


t  kllxlhm  cos(  kllxlhm  hl)°os(  kiix2hm  h2^  “ 

knx2hm  sin<  kIIxlhm  hi>sin<  kllx2hm  h2>  1  2t°-5h3  ' 
0.25sin  (2  kIIx3hm  h3)/  kIIx3hm  ]  +  [sin(  knxlhm  h^]2 

[  0.5h2  +  0.25sin  (2  knx2hm  h2)/kIIx2hm  3  + 

f  kIIxlhm  cos(  kIIxlhm  ^1^  /,kIIx2hm  J2 
[  0.5h2  -  0.25sin  (2  kIIx2hm  h 2)/  knx2hm  ]  + 

sin(  kjixihm  *h)  kllxihm  cos(  kllxlhm  ^l)  /  kllx2hm 

i  sin<  knx2hm  m2  ym 


( C-15) 


i 


■  * .*».»  ut  1. 1  ■»«  i 


Ntnim  ^  ^  ^  0.5hj  +  0.25sin  (2  kjixierrl  hj)/  kllxlem  ]  /  Ej-t 

f^cosC  kIIx3em  h3)]'2  [  EjCOs(  kIIXlem  hj)cos(  kIIx2em  h2)  - 

^^Ilxlem  sm(  kIIxlem  hi)sin(  kiix2em  h2^  t  ei  kllx2em  J  2/  h 
[0.5h3  +  0.25sin  (2  kjjx3em  h3)/  kjIx3em  ]  + 

f  ‘kIIxlem  ^sin(  knxlem  ^i)  /  kllx2em  ^  !  ^2 
[  0.5h2  -  0.25sin  (2  k„x2em  h2)/  kIIx2em  ]  + 

[  EjCos(  kIIXlhmh1)/e2]2/E2 
[  0.5h2  +  0.25sin  (2  kjjx2hm  h2)/  kjlx2jim  ]  + 

J  *knxlem  ^sin(  kjjxiern  h|)  /  kjIx2em  h2]  [  EjCOs(  kj^jj^  hj) 


2/r  1-1/2 


[  sin(kIIx2hmh2)]^/e2  } 


( C-16 ) 


The  othercoefficients  such  as  Ntjj2m<  Nrjj2m  can  be  obtained  from  Eqs.  (C-6) 
to  (C-13). 
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APPENDIX  D.  DERIVATION  OF  THE  SCATTERING  MARTTX 
FOR  THE  DOUBLE-STEP  JUNCTION 


In  reference  to  Fig.  II-3  and  Eqs  ( H-3)  to  ( 13-24),  the  tangential  fields  in 
Region  I  can  be  written  as 

Ely  =  -  d<j>j  /dz  +  (1  /jcoeQ)  dfyj  /  dxdy 
M  N-i 

=  -  S  X  ^izhmn^'^Imn  +  ®Imn)  ^Imn ' 
m=l  n=0 
M-l  N 

(l/jcoe0)ZZ  N Iemn (mn Jt^/a j b j )  sin(nmx/ai)  cosCrnty/bj)  (Cjmn  +  DImn) 

m=0  n=l 

(D-l) 

EIx  =  (1  /  jcoe0)  (  d2  /  dx2+  kr2  )  \j/j 

M-l  N 

=  (l/jcoe0)  Z  Z  [  kQ2-  (mn  /  ai)2]  (Cj^-t-  Djj^)  Qjmn 

m=0  n=l 

(D-2) 

Hjy  =  dy{/dz  +  (1  /j(0(Io)  92<J>j  /  dxdy 
M-l  N 

=  Z  Z  klzemn^Imn  +  ^Imn^  Qlmn  " 

m=0  n=I 

M  N-l 

(1  /  jo)Hc)  Z  Z  N^^mnTAa^)  cosCmroc/aj)  sinCmty/b!)  (AImn+BImn) 
m=l  n=0 

(  D-3  ) 

HIx  =  (1  /  jcoja.0)  ( d2  /  9x2+  kj2 ) 

M  N-l 

=  (l/j(0iio)  Z  Z  [  kQ2-  (mn  /  ai)2]  (AImn+  BImn)  PIinn 
m=l n=0 

( D-4  ) 


where 
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Almn  =  &  TmnexP  (~kl7hmnz)  «  ®Imn  “  ^  TmncxP  (kT7.hmnz) 


^Imn  =  cImnexP  ("kIzemnz)  *  ^Tmn  -  d  TmnexP  (kLzemnz) 


( 

In  the  Region  II  the  tangential  fields  are  given  by 

Ejjy  =  -  5<f>n/3z  +  (1  /jcoe0£j)  d\{  /  dxdy 
p  Q-l 

=  -  L  Z  kHzhpq('^IIpq  +  ®IIpq)  EIIpq  +  O/j^o^vIIq' 

p=l<t=0 
p  Q 

Z  Z(1  /€j  (X))  (qn/b)cos[q7t/b(y-y1)]  (-Cnpq  +  Dnpq)  dT^  dx 
p=0q=l 

( 

eHx  *  C1  /j“e0ej)  ( d2  /  9x2+  knj2 )  yn 
P-1Q 

=  knj2-  kHxjep2^  ^Hpq+  ^Ilpq)  Qllpq  !  ej(x) 

p=Oq=l 


Hny  =  dvn  /  3z  +  (1  /  j(OU0)  32<Pn  /  dxdy 
P-1Q 

=  Z  Z  kHzepq(~^IIpq  +  ^Ilpq)  Qllpq ' 

p=0q=l 

P  Q 

(1  /  jcon^  Z Z  NIIsq(q7t/b)sin[q7t/b(y-y 2)]  (AIIpq+  Bnpq)3Rp/  dx 
p=l  q=0 

Hnx  =  0  /  < 92 ' 9x2+  kn2 )  ^ 

p  Q 

=  (1/jcop.Q)  Z  Z  ( knj2-  kIIxjhp2)  ^Hpq+  EIIpq)  EIIpq 


where 


^Ilpq  =  a  IIpqexP  (‘kllzhpq2)  i  ®Hpq  “  ^  DpqexP  ^Ilzhpq2) 


^Ilpq  _  cnpq^P  ("kllzepq2)  *  ^Ilpq  “  d  npqexP  (^Ilzepq2) 


Match  the  tangedal  fielde  at  z  =  0  and  in  area  A2  leads  to 


(D-10) 


M-l  N 


I  I[kG2-  (mn/  ai)2]  (Clmn+  dImn) 

m=On=l 
P  Q 

=  Z  Z[  (  qit/b)2-  kjlzjep2]  (cIIpq+  ^Ilpq)  Qllpq/£j(x) 
p=Oq=l 


(D-ll) 


M  N-l 

^  ^  ^Izhmn^'aImn  +  ^Imn^  ^Imn  + 
m=l  n=0 

M  N 

(l/jaJEoJZ  I  NjenjnCmnAjbj)  suKnmx/aj)  cosCnrcy/b!)  (cTtT(n  +  dj^) 

m=0 n=l 
P  Q-l 

=  Z  Z  kllzhpq('allpq  +  ^Ilpq)  ^Ilpq '  (1/jweQ^  ^vIIq' 
p=lq=0 
P-1  Q 

Z  Z  (qrc/b)cos[q7t/b(y-y j )]  (-cnpq  +  dnpq)  dTJ  dx  /  Ej(x) 
p=Oq=l 

(D-12) 


M  N-l 


Z  Z  [  ko  (m7C  /  ai)  1  (aImn+  ^Imn^  ^Imn 


m=l  n=0 


-  Z  Z  [  (  qtt/b)  -  kjj2jep  ]  (ajjpq+  bjjpq)  Pjjpq 

p=l  £p=0 


(D-13) 
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2  Z  ^IIzepq('cImn  +  dImn'  Qlmn " 
m=0  n=l 

MN-1 

(1  / jflHijj)  1 1  NjhnjnCmnAibj)  cosCnmx/ai)  sinCnTcy/bj)  (ajmn+bimn) 
m=l n=0 

P-l  Q 

=  £  Z  ^Hzepq^IIpq  +  dIIpq)  Qllpq  ‘ 

psOqsl 

P  Q 

(1  /  jcoUq)  Z  Z  Nsnq(q?^b)sin[qn/b(y-y  j )]  (anpq+  bnpq^V  dx 

Dal  0=0 

(  D-14  ) 

After  applying  the  orthogonality  of  the  eigenfimetins,  following  equations  are  obtained 


M  N-i  M  N_1 

Z  Z  Mai(mji,p,q)aIinn-Vbn(p,q)bnpq=  -Z  Z  MaI(m,n,p,q)  bj^  +  Vbn(p,q)  anpq 

maljiaO  m-l.i*«0 

(D-15  ) 

P-l  Q  P-l  Q 

-Z  Z  Mdn(m,n,p>q)dnpq+Vci(m,n)cim  =  Z  ZM(jjj(m,n,p,q)  cjjpq  -  VcI(m,n)  djmn 

paO,<fi  p=0.q=1 

(D-16) 

M  N-l  M-l  N 

Z  Z  MbI(m,n,p,q)aImn+Z  Z  McI(m,n,p,q)  -VaII  (p,q)  bIIpq  +  Vdn(p,q)  dnpq 

m=lji=0  m«0,n=l 

MN-1  M-l  N 

=  -Z  Z  MbI(m,n,p.q)bImn+Z  Z  Mcl(P*n.P*£l)dImn+van(P'tl)aIIpq+  VdIl(P>cO  cIIpq 

maljlaO  fn*0Jl*l 

( D-17) 


PQ-l  P-l  Q 

Z  Z  Mbn(p,q^n,n)  bnp+I  Z  MdjCp.q.m.n)  dflpqVai  (mfn)aIm-VdI(m,n)  cj^ 
p*l,q-0  P=0.q=l 

PQ-l  P-l  Q 

=Z  Z  Mbu(p,q^n,n)  ajjp-Z  Z  Mcji(p,q,m,n)  cjjpq+Vjj  (m,n)  bimn+VdJ(m,n)  dj^ 


where 


yj+bX!** 

MjjCm.n.p.q)  =  [  l^2-  (mx  /  aj)2]  \  dy  \  dx  Pimn(x,y)  Pnpq(x*y) 

yi  *1 

y^xj+x 

Mdn(m,n,p,q)  =  [  ( q*/b)2-  kj^^2]  f  dy  J  dx  Qtan(x,y)  Qnpq(x*y)  /  ej« 

yi  *1 

Mbj(m,n,p,q)  yi^*i+* 

=  Njemn )  ^y  J  ^  Qlmn^x,y^  Qllpq^x*y)  I  ejM 

yi  *1 

yi+*>  *i+* 

McI(m,n,p,q)  =  J  dy  J  dx  CWx.y)  Qiipq(x,y)  /6j(x) 

yi  *1 

y,+b  Xi+x 

Mbn(m,n,p,q)  =  kjj^pq  J  dy  I  dx  Qinm(x.y)  Qnpq(x>y) 

yi  xi 

Mcn(m,n,p,q)  yi+**i« 

=  Nvjjqqx/  0weob)  I  dy  J  dx  cos[qx/b(y-yj)  Pimn(x»y)^Tnp(x)  /  ^x  /  ej  M 
yi  *i 

vbn(p>q)  =  ( q^)2-  k2nzepq 
VcI(m,n)  =  kQ2-  (mx  /  aj)2 
=  (  qx/b)2-  k2nZCpq 
van(P-q)  yi+b 

=  NjUqqx/  (jd))i0b)  /  dy  /  dx  sin[qx/b(y-y  j )  Qiipq(x.y)^Riip(x>)  /  ^x  /  ej  (x) 
yi  *i 


vbn(P*q) 88  ^Ilzepq 
val^m,n^  “  ^Ilzhtnn 


VdjCnwO  =  mnjAcajb!  jcoe0 ) 

Equations  ( D-15 )  to  ( D-18 )  may  be  written  in  matrix  form 


Maj  0  -Vbn  0  ' 

AImn' 

0  vcl  o  -Mdn 

Mmn 

MblMcI  'vanvdn 

BDpq 

val"vdl  Mbn 

f-Ma!  0  Vbn  0 

'  BImn 

1 

O 

> 

• 

e 

BImn 

-Mbj  Mcj 

AIIpq 

Val  Vd j  Mbn  -McII 

Cnpq , 

( D-19 


where  is  a  matrix  of  dimension  of  PxQxMxN  and  is  given  by : 


MjjU, 0,1,0),  M^l.0,1,1) ...  MaI(l,0,M,N-2),  M^lAM^-l)" 
Mjjfl, 1,1,0),  Mjjd.l.U) ...  MaI(l,l,M,N-2),  Mjjd.l.M.N-l) 


= 


MaI(2, 0,1,0),  MaI(2, 0,1,1) ...  MaI(2,0,M,N-2),  M^.O.M.N-l) 
Maj(2, 1,1,0),  MaI(2,l,l,l) ...  MaI(2,l,M,N-2),  MaI(2,l,M,N-l) 


MjjCP.Q-2,1,0),  MaI(p,Q-2,l,l) ...  MajC..),  MaI(P,Q-2,M,N-l) 

.  Mjj(P,Q- 1,1,0),  MjjCp.Q- 1,1,1) ...  Mjjf...),  MaI(P,Q-l,M,N-l)  , 

The  dimensions  of  matrix  Mjjjj  Mbj  Mcj  >  M^jj  and  Mcu  are  PxQxPxQ, 


MxNxMxN,  MxNxMxN,  MxNxPxQ,  and  MxNxPxQ,  respectively.  The  Vbjj  is  a 
diagonal  matrix  of  PxQxPxQ.  It  is  represented  by 


0 

0 


vn(0,0)  0  0  .  .  . 

o  Vt,u(l,l)  0  .  .  . 

0  Vbn<2,2) 


0 

0  0  0 


0 

0 

0 


0 

Vbn(Q-2,Q-2)  0 

•  Vbn(Q-l,Q-l) 


Other  diagonal  matrices  such  as  Vc^  Vbjj  ...etc.  may  be  obtained  similarly. 
In  Eg.  (D-19),  Ajjjjjj ,  Djipq  •  •••  are  the  coefficient  vectors  of  dimensions  of 
MxN,  Px  Q,...  and  so  on.  For  instance,  Ajmn  is  a  vector  of  dimension  of  MxN  and 
is  givenby 

AImn  =  £aI10,aIll,  •  aI20,  aI21,  •>  aIm(n-2)  >  aIm(n-l)] T 
where  T  indicates  a  transpose  operation. 


Ok. 


APPENDIX  E.  DERIVATION  OF  SCATTERING  MATRIX  FOR 


THE  CASCADED  DISCONTINUITIES 


From  Fig.  II- 10,  the  following  matrix  equations  can  be  derived: 


R(A)j  fSfin  Sflujfiw 

H  '  SfUl  Sf]22  I®]  <  E-l ) 

>)  |  (o  S,-)  I<2) 

I<3>  J  st  0  1(3).  ( E-2 ) 

'  /  / 

R®'  [son  s„i2j  f  '<3) 

R<«j  ‘  Sfl21  SH22  M  (E-3) 

/ 

Equation  ( E-2 )  is  substituted  into  ( E-l )  to  obtain 

R<a>  =  SnilI<A>  +  Sf 1 12  S,  R(»  (E-4) 

R®  -  SflIlI<A>  +  Sfl22  S,  R®  (E-5) 

Next,  ( E-2 )  is  substituted  into  ( E-3 )  to  get 

R(3)  =  S^n  St  RC2)  +  Snn  KQ  ( E-6 ) 

R(C)  =  Sf221  St  R<2>  +  Sf222 1(Q  ( E-7  ) 


'  *  a.  ¥  *  w  ■  J*  a  "  a."  » 


’  /  V  V 

JuAAAtuTjiTj’j  v  •-»  v  fj>  ’J  *_*  r. 


.'O 


Equations  ( E-5)  and  ( E-6 )  are  used  to  isolate  R®  and  Rw  : 

R©  -  U,  Sn21  I<AH  U,  Sn22  S,  S(2,2 1(C'  (  E-8 

R<»  =  U2  Sf2U  St  Stl2I  I<A>  +  U2  Sf2,2  I<A>  (  E-9 

where 

ui  =(i-  sn22stsnilsty'  (B-l 

U2  =  (I -  Sf2n  st  Sfi21  st y1  ( E-l 

Where  I  is  the  identity  matrix.  Finally,  equations  ( E-8)  and  ( E-9)  are 
substituted  into  equations  ( E-4 )  and  ( E-7) 

R<a>  =  [Snn+Sni2S,U2SGlIS,  Sfl21]  1(A) 

+  Sfll2S,U2  SqijI®  ( E-l 

r<q  -  sC21  S,  U,  Sn21I<A>  + 

l  Sf222  •*•  Sf221  StUl  Sfl22  St  Sf2l2l  I<C)  (  E-l 

The  submatrixes  S^,  SAC>  SCA,  and  Scc  are  easily  idetified  from  the  above 
equations. 
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